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In this paper we study the possibility of testing Charge-Parity-Time Reversal (CPT) symmetry 
with cosmic microwave background (CMB) experiments. We consider two kinds of Chern-Simons 
(CS) term, electromagnetic CS term and gravitational CS term, and study their effects on the CMB 
polarization power spectra in detail. By combining current CMB polarization measurements, the 
seven-year WMAP, BOOMERanG 2003 and BICEP observations, we obtain a tight constraint on 
the rotation angle Aa — —2.28 ± 1.02 deg (ler), indicating a 2.2 a detection of the CPT violation. 
CN| . Here, we particularly take the systematic errors of CMB measurements into account. After adding 

the QUaD polarization data, the constraint becomes —1.34 < Aa < 0.82 deg at 95% confidence level. 
When comparing with the effect of electromagnetic CS term, the gravitational CS term could only 
\ generate TB and EB power spectra with much smaller amplitude. Therefore, the induced parameter 

(— i ■ e can not be constrained from the current polarization data. Furthermore, we study the capabilities 

" of future CMB measurements, Planck and CMBPol, on the constraints of Aa and e. We find that 

the constraint of Aa can be significantly improved by a factor of 15. Therefore, if this rotation angle 
effect can not be taken into account properly, the constraints of cosmological parameters will be 
■ biased obviously. For the gravitational CS term, the future Planck data still can not constrain e very 

well, if the primordial tensor perturbations are small, r < 0.1. We need the more accurate CMBPol 
experiment to give better constraint on e. 



PACS numbers: 98.80.Cq, 98.80.Es, 11.30.Cp, 11.30.Er 
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Q-[ I. INTRODUCTION 

In the standard model of particle physics, Charge-Parity-Time Reversal (CPT) symmetry is a fundamental sym- 
c/3 ■ metry. Probing its violation is an important way to search for the new physics beyond the standard model. Up to 
now, CPT symmetry has passed a number of high-precision experimental tests and no definite signal of its violation 
has been observed in the laboratory. So, the present CPT violating effects, if they exist, should be very small to be 
amenable to the experimental limits. 

However, the CPT symmetry could be dynamically violated in the expanding universe. In the literature [H-Q, the 
r*"** ■ cosmological CPT violation mechanism has an interesting feature, which is that the CPT violating effects at present 
are too small to be detected by the laboratory experiments, but large enough in the early universe to account for the 
generation of matter-antimatter asymmetry. More importantly, these types of CPT violating effects could be accu- 
mulated to be observable in the cosmological experiments (j-Q ■ With the accumulation of high-quality observational 
data, especially those from the cosmic microwave background (CMB) experiments, cosmological observation becomes 
' a powerful way to test the CPT symmetry. 

In this paper we firstly consider the cosmological CPT violation in the photon sector where the electrodynamics 
is modified by a Chern-Simons term £ ecs ~ p ll A v F ,J,v . Here, p^ is an external vector and F^ v = (l/2)e A1!ypo 'F pa is 
the dual of the electromagnetic tensor. This term violates the Lorentz and CPT symmetries if p^ is treated as an 
external field. One of the physical consequences of this electromagnetic Chern-Simons (eCS) term is the rotation of the 
$_i ' polarization direction of electromagnetic waves propagating over large distances Q , which is known as "cosmological 
birefringence" . This rotation angle Aa can be obtained by observing polarized radiation from distant sources such as 
radio galaxies, quasars, and CMB. For the standard theory of CMB, the TB and EB cross-correlation power spectra 
vanish. In the presence of the non-zero rotation angle Aa, one expects to observe non-zero TB and EB power spectra, 
even if they are zero at the last scattering surface. Denoting the rotated quantities with a prime, one gets [H, Q 
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while the CMB temperature power spectrum remains unchanged. Hence we can use the CMB polarization measure- 
ments to test the Lorentz and CPT symmetries in this eCS model (see Refs. H4H1 , and references within) . With 
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homogeneous and isotropic rotation angle, Ref.[15| used the newly released CMB observations and found that a non- 
zero rotation angle Aa = —2.2 ± 0.8 deg (1 a) is favored by the CMB polarization data from the seven-year WMAP 
(WMAP7) [3, BOOMERanG 2003 (B03) and BICEP [H observations. 

Besides the eCS term, the gravitational Chern-Simon s (g CS) term could also generate the non-zero CMB TB and 
EB cross-correlation power spectra. Following Refs. |2lTl23l |. here we consider a Lorentz and CPT violating term in 
the gravity sector £ gcs ~ ^ v p° ' R'p R^ pa , where R'su.u IS tnc Riemann tensor. This gCS term does not affect the 
evolution of background and scalar perturbations. Therefore, the effect of gCS term only appears in the evolution of 
tensor perturbations, if we neglect the vector perturbations. As we know, the gravitational wave has two independent 
polarized components denoted by + and x. We usually use the right- and left-handed circular polarized components: 

h R = -j=(h+-ih*) , h L = ^=(h + +i hX) . (2) 

The power spectra for different handedness are defined as: (h s *(ki)h s (It?)) = P^ s <5 3 (ki — k2)<5 SS ', where superscripts 
s and s' stand for the circular polarized state, s, s' = R, L. Due to the presence of this gCS term, the produced power 
spectra P R and Pj^ are not equal. Similarly with Refs.|2fT[23|. we use a parameter e to characterize this discrepancy: 

Ph = \W - e) , P^ = \p h {l + e), 

Pf + Pl=P h , P* - P£ = -eP h . (3) 

We can see that e = — 1, 0, 1 denote purely right-handed polarized, unpolarized and purely left-handed polarized 
gravitational wave, respectively. In this paper we only consider the scale-independent e simply. 
The CMB power spectra generated by the tensor perturbations are given by: 

Cf x ' = (At:) 2 J k 2 dk[P R (k)+P^(k)]Af(k)Af' (k) = (An) 2 J k 2 dkP h (k)Af(k)Af'(k), 

Cj Y ' = (An) 2 f k 2 dk[Pl(k) - Pk(k)]Aj(k)Aj'(k) = -(An) 2 , f k 2 dkP h (k)Aj (k)Aj' (k) , (4) 



where XX' and YY' denote TT, TE, EE, BB and TB, EB, respectively, and Af(k) are the transfer functions. These 
equations show that the TT, TE, EE and BB power spectra only depend on the sum of the primordial power spectra 
of gravitational waves, while the TB and EB power spectra rely on the difference between the power spectra of right- 
and left-handed polarized components —ePh- Therefore, in the presence of the gCS term, the non-zero TB and EB 
correlations will be generated and other four power spectra are unchanged. 

In this paper we will study the effects of eCS and gCS terms on the CMB power spectra in detail and perform a 
global analysis on them by using the latest CMB polarization measurements, as well as the future simulated CMB 
data. The structure of the paper is as follows: in section [TT] we describe the current and future simulated datasets we 
use. Section \TU\ contains our main results from the current observations and future measurements, while section UVl is 
dedicated to the conclusions and discussion. 



II. CMB DATASETS 



A. Current Datasets 

In our calculations we mainly use the full data of WMAP7 temperature and polarization power spectra (HI- The 
WMAP7 polarization data are composed of TE/TB/EE/BB/EB power spectra on large scales (2 < t < 23) and 
TE/TB power spectra on small scales (24 < I < 800), while the WMAP7 temperature data are only used to set the 
underlying cosmology. For the systematic error, the WMAP instrument can measure the polarization angle to within 
±1.5 deg of the design orientation (2~Ij ]. In the computation we use the routines for computing the likelihood supplied 
by the WMAP team. Besides the WMAP7 information, we also use some small-scale CMB observations. 

The BOOMERanG dated January 2003 Antarctic flight measures the small-scale CMB polarization power 
spectra in the range of 150 < I < 1000. Recently, the BOOMERanG collaboration re-analyzed the CMB pow er 
spectra and took into account the effect of systematic errors rotating the polarization angle by —0.9 ± 0.7 deg |25| . 

Recently, the Background Ima ging of Cosmic Extragalactic Polarization (BICEP) [20j and QU Extragalactic Sur- 
vey Telescope at DASI (QUaD) [26[ collaborations released their high precision data of the CMB temperature and 
polarization including the TB and EB power spectra. These two experiments, locating at the South Pole, are the 
bolomctric polarimcters designed to capture the CMB information at two different frequency bands of 100GHz and 



3 



TABLE I: Assumed experimental specifications. We use the CMB power spectra only at I < 2500. The noise parameters At 
and Ap are given in units of /iK-arcmin. 



Experiment 


/sky 


t-max 


(GHz) 


#FWHM 


A T 


Ap 


PLANCK 


0.65 


2500 


100 


9.5' 


6.8 


10.9 








143 


7.1' 


6.0 


11.4 








217 


5.0' 


13.1 


26.7 


CMBPol 


0.65 


2500 


217 


3.0' 


1.0 


1.4 



150GHz, and on small scales - the released first two-year BICEP data are in the range of 21 < I < 335 [20j ; whereas 
the QUaD team measures the polarization spectra at 164 < I < 2026, based on an analysis of the observation in the 
second and third season [H, [l3|. They also provide the systematic errors of measuring the polarization angle, ±0.7 
deg and ±0.5 deg, for BICEP and QUaD observations, respectively. 



B. Future Datasets 



Current CMB measurements are still not accurate enough to verify the possible CPT violation. We follow the 
method given in Refs.0, [TtJ and simulate the CMB power spectra with the assumed experimental specifications of 
the future Planck [53] and CMBPol [28| measurements. We choose the best-fit model from the WMAP7 data [lH as 
the fiducial model. 

In Table U we list the assumed experimental specifications of the future Planck and CMBPol measurements. The 
likelihood function is C oc cxp(— x^ ff /2) and 

x t s = J2(2£ + l)/ sky (^ + In jgj + 3j , (5) 

where / s k y denotes the observed fraction of the sky in the real experiments, A is defined as: 

a = offerer - (crr)+cj E (crcr -cj E cn 

, riTB i r<TE r<EB r<TBr<EE\ , r^TE i f^TB f^EB f^TE f^BB\ 
-t- O f Ly ( — Ly ( L> e ) -t- O f L> e — L> e O f ) 

+ C EE (Cj T C BB - {CJ B f) + C EB {CJ E Cj B - Cj T C EB ) 

i nTB I r<TE r<EB FiEE FtTB \ | r<EB i r<TE riTB r<TTr,EB\ 
-t- (O^ — L> £ ) -t- £ \Kj t O f — £ O f ) 

+ C BB {Cj T C EE ~{Cj E f), (6) 
and |C| and |C| denote the determinants of the theoretical and observed data covariance matrices respectively, 

i/=fi _ rtTT riEE r<BB , 0/ ^,TE /=tTB r<EB rtTT(r<EB\2 

l°l — °£ W W ' & W W — u i \ye. ) 
-C EE {Cj B f-C BB {CJ E )\ 

\r<\ _ r<TTriEEr<BB , o/^TErtTBriEB ATT/pEB\2 

M — W W W + L( -"i W — W V^e. ) 

-C EE {Cj B Y~C BB (Cj E f . (7) 

The likelihood has been normalized with respect to the maximum likelihood XcS = 0, where Cf" Y = Cf~ Y . 



III. NUMERICAL RESULTS 



In our study we make a global analysis to the CMB data with the public available Markov Chain Monte Carlo 
package CosmoMC [29j . which has been modified to compute the non-zero TB and EB power spectra discussed above. 
We assume the purely adiabatic initial conditions and impose the flatness condition motivated by inflation. Our basic 
parameter space is: P = (ut>, w c , J7a, t, n s , A Sl r), where wj, = f2f,/i 2 and uj c = £l c h 2 are the physical baryon and cold 
dark matter densities relative to the critical density, fl\ is the dark energy density relative to the critical density, r 
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TABLE II: Constraints on the rotation angle Aa (68% C.L.) from various CMB data combinations. 



CMB Data sets 


without systematics 


with systematics 


WMAP7 


-1.06 ±1.39 deg 


-1.14 ±2.05 deg 




— o.y / ± 4.U5 deg 


A ao J_ A 1 a . ] .. . 

—4. bo ± 4.1b deg 


BICEP 


-2.59 ± 0.99 deg 


-2.57 ± 1.21 deg 


QUaD 


0.60 ± 0.40 deg 


0.59 ± 0.64 deg 


WMAP7+B03+BICEP 


-2.23 ±0.83 deg 


-2.28 ± 1.02 deg 


WMAP7+BICEP+QUaD 


0.07 ± 0.36 deg 


-0.22 ± 0.54 deg 


WMAP7+B03+BICEP+QUaD 


0.04 ± 0.34 deg 


-0.26 ± 0.54 deg 



is the optical depth to re-ionization, A s and n s characterize the primordial scalar power spectrum, r is the tensor 
to scalar ratio of the primordial spectrum. For the pivot of the primordial spectrum we set k s o = 0.002 Mpc -1 . 
Furthermore, in our analysis we include the CMB lensing effect, which also produces B modes from E modes [30| . 
when we calculate the theoretical CMB power spectra. 

A. Electromagnetic Chern-Simons Term 

Firstly, we consider the constraint on the rotation angle Aa, induced by the eCS term, from the current CMB 
measurements. As we know, this rotation angle is accumulated along the journey of CMB photons, and the constraints 
on the rotation angle depends on the multipoles I [l(|. Refs.fltlllii found that the rotation angle is mainly constrained 
from the high-^ polarization data, and the polarization data at low multipoles do not affect the result significantly. 
Therefore, in our analysis, we assume a constant rotation angle Aa at all multipoles. Further, we also impose a 
conservative flat prior on Aa as, — tt/2 < Aa < ir/2. 

In our previous works, we only presented the statistical errors of rotation angle and did not consider the possible 
systematic errors of CMB measurements. Inspired by Ref . [25| . in this paper we consider two rotation angles, Aa and /?, 
in order to take into account the real rotation signal and a systematic error for each CMB polarization measurement. 
Therefore, in our analyses we have five free parameters in this eCS model: the rotation angle signal Aa and four 
systematic errors, /3wmap7, Pbo3, /3bicep, /?QUaD, for four CMB observations, respectively. And we impose priors on 
these four systematic errors: 

/?wmap7 = 0.0 ± 1.5 deg , /3 B03 = -0.9 ± 0.7 deg , 
/3bicep = 0.0 ± 0.7 deg , /3 QUaD = 0.0 ± 0.5 deg , (8) 

and marginalize over them to constrain the rotation angle. In Table [TT] we present current constraints on Aa from the 
WMAP7, B03, BICEP and QUaD CMB polarization power spectra with and without CMB systematic errors. 

When only using WMAP7 power spectra data at all multipoles £, without the systematic effect (/?wmap7 = 0), 
we obtain the constraint on the rotation angle: Aa = —1.06 ± 1.39 deg at 68% confidence level, which is quite 
consistent with that obtained from the WMAP team [l4j and is a significant improvement over the WMAP3 [| and 
WMAP5 [l(J UK results. After adding the prior of /3wmap7 into the x 2 calculation, the constraint becomes weaker: 
Aa = -1.14±2.05 deg (68% C.L.) and -5.24 < Aa < 2.96 deg (95% C.L.). Similarly, we revisit the constraint on Aa 
from the B03 polarization data and obtain the constraints at 68% C.L.: Aa = — 5.97±4.05 deg and Aa = — 4.63±4.16 
deg, without and with the CMB systematic effect, respectively. The results are in good agreement with the previous 
results of Refs.0, [25|. These two CMB measurements show a weak indication for a non-zero rotation angle about 1 a 
confidence level. 

Recently, BICEP and QUaD collaborations also released their high precision CMB polarization data. In our previous 
work [TH, we found that the BICEP data alone could give very tight constraint on the rotation angle: Aa = — 2.59±0.99 
deg (68% C.L.), when fixing /3bicep = 0, which means that BICEP data alone favors a non-zero rotation angle at 
about 2.5(7 confidence level. When we include the systematic effect of BICEP measurement, the constraints at 68% 
and 95% C.L. are: 

Aa = -2.57 ± 1.21 deg (68% C.L.) , - 4.99 < Aa < -0.15 deg (95% C.L.) , (9) 

which still gives a more than 2 a detection of a non-vanishing rotation angle. Furthermore, when WMAP7 and B03 
data are added to the BICEP sample, the constraint on Aa gets tightened: 



Aa = 



-2.28 ± 1.02 deg (68% C.L.) , 



- 4.32 < Aa < -0.24 deg (95% C.L.) , 



(10) 
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FIG. 1: Left panel: The one-dimensional posterior distributions of the tensor to scalar ratio r derived from the 
WMAP7+B03+BICEP data combinations including (black solid line) and not including (red dashed line) the rotation an- 
gle. Right panel: Two-dimensional 1,2 a contours in the Aa — r plane. 



which implies Acv ^ at 2.2 a confidence level, even considering systematic effects of these three CMB measurements. 

Finally, we constrain the rotation angle from the QUaD polarization data. Similarly with previous results of QUaD 
collaboration, we use the QUaD data alone and obtain the constraint at 68% confidence level without and with 
systematic effect: Aa = 0.60 ± 0.40 deg and Aa = 0.59 ± 0.64 deg, respectively. When comparing with the result of 
WMAP7+B03+BICEP [Eq.dTUJ], there is a - 2 a tension between QUaD and WMAP7+B03+BICEP observations, no 
matter whether we include the systematic effects of these CMB measurements. As we discussed before [Hj], this tension 
should be taken care of in the further investigation. By combining WMAP7, BICEP and QUaD data and including 
their systematic effects, we obtain the constraint: Aa = — 0.22 ± 0.54 deg (68% C.L.) and —1.30 < Aa < 0.86 deg 
(95% C.L.), which is consistent with the result of Ref . (l4| . When adding the B03 data into the analysis, the constraint 
on the rotation angle does not change significantly, due to large error bars of B03 polarization data. 

Moreover, the sources of the CMB polarization, especially for the B-mode, are not unique. For example, the B- 
mode can be generated by the cosmological birefringence as mentioned above; it might be converted from E-mode by 
cosmic shear [30j; it could be the signature of the gravitational waves; or it can even be produced by the instrumental 
systematics. Therefore, one should bear in mind that the rotation angle might be degenerate with other cosmological 
parameters or nuisance parameters when fitted to the polarization data. Therefore, in order to distinguish these 
effects and obtain the clean information of the primordial tensor B-mode, the rotation angle has to be constrained, 
and the measurements of TB and EB power spectra are really necessary. Otherwise, the constraints of cosmological 
parameters might be biased, if this cosmological birefringence effect can not be taken into account properly. In the 
left panel of FigJTJwe show the one-dimensional posterior distributions of the tensor to scalar ratio r derived from the 
WMAP7+B03+BICEP data combinations. When including the rotation angle, the 95% upper limit of r is: r < 0.26, 
while we obtain r < 0.29 with Aa = 0. Since the data combination WMAP7+B03+BICEP favors a negative rotation 
angle which contributes the BB power spectrum, the BB power spectrum from primordial tensor perturbations, as 
well as the upper limit of r, will be slightly suppressed, consequently. We also show the two-dimensional contour 
between the tensor and the rotation angle in the right panel of Fig[TJ Due to the large error bars of current CMB 
polarization data, the effect of non-zero rotation angle is neglectable on the constraints of cosmological parameters. 
However, this effect might be important in the future forecasts. 

Therefore, we simulate the future CMB power spectra with Planck to investigate the effect of non-zero rotation 
angle on the constraints of cosmological parameters. The fiducial model we choose is the best-fit WMAP7 model [Til ]: 
fl b h 2 = 0.023, fl c h 2 = 0.108, fl A = 0.752, r = 0.087, n s = 0.98, and log [IQ W A S ] = 3.135 at k s0 = 0.002 Mpc" 1 . 
For the tensor to scalar ratio r, we have two fiducial models: r = and r = 0.3. For each model, we consider three 
fiducial values of the rotation angle with Aa = 0, —2, —5 deg. Here, we neglect the systematic error of future CMB 
measurement and the CMB lensing effect. In Table IIIII we list the constraints on the cosmological parameters with 
different fiducial models. 

When we include the rotation angle in the calculations, as you can see in Table IIIII the fiducial values of input 
parameters are always recovered. The future Planck measurement could shrink the standard deviation of rotation 
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TABLE III: Constraints on the cosmological parameters (68% C.L.) from future CMB mock data. 



Parameters 


100 Q b h 2 


10 tt c h 2 


10 SI a 10 T 10 n s 


loefl0 10 /Ul 


r(95% C.L.) 


/\rv fHpo"^ 


Ay 2 


Fiducial 


2.300 


1.080 


7.520 0.870 9.800 


3.135 


See below 


See below 


o 










PLANCK, r = 0, Aa free 










Aa — 


no 
U 


2 300 + 017 


1 079 + 015 


7 519 + 077 868 + 057 9 804 + 045 


3 1 33 + 01 9 


< 0.033 


—0 001 + 058 


o 


A ... 


—2 


2.301 ± 0.017 


1.079 ± 0.015 


7.521 ± 0.076 0.870 ± 0.056 9.805 ± 0.043 


3.133 ± 0.019 


< 0.035 


—2.001 ± 0.061 





\ ■ = 


—5° 


2.301 ± 0.017 


1.079 ± 0.015 


7.520 ± 0.078 0.867 ± 0.057 9.805 ± 0.045 


3.133 ± 0.019 


< 0.037 


—5.000 ± 0.062 


o 










PLANCK, r = 0, Aa = 










A . . 

Aa = 


0° 


2 300 + 017 


1 079 + 01 5 


7 518 + 078 870 + 057 9 806 + 044 


3 133 + 019 


< 0.032 




o 


A ... 


—2° 


2.304 ± 0.017 


1.082 ± 0.015 


7.506 ± 0.080 0.868 ± 0.056 9.804 ± 0.045 


3.134 ± 0.019 


< 0.048 




6.98 


\ — 


-5° 


2.327 ± 0.017 


1.085 ± 0.015 


7.499 ± 0.077 0.858 ± 0.056 9.836 ± 0.044 


3.119 ± 0.019 


0.081 ± 0.058 




261.68 










PLANCK r = 0.3 Aa free 










A« = 


0° 


2.300 + 0.017 


1.080 + 0.015 


7.512 + 0.078 0.870 + 0.056 9.801 + 0.045 


3.135 + 0.020 


0.305 + 0.100 


0.003 ± 0.060 





Aa = 


-2° 


2.300 + 0.016 


1.080 + 0.015 


7.513 + 0.077 0.871 + 0.056 9.802 + 0.044 


3.135 + 0.020 


0.308 + 0.100 


-2.003 + 0.059 





Aa = 


-5° 


2.300 + 0.017 


1.080 + 0.015 


7.513 + 0.076 0.873 + 0.057 9.800 + 0.044 


3.136 + 0.020 


0.303 + 0.102 


-4.998 ± 0.059 













PLANCK, r = 0.3, Aa = 










Aa = 


0° 


2.300 + 0.017 


1.080 + 0.015 


7.515 + 0.078 0.872 + 0.057 9.800 + 0.045 


3.135 + 0.020 


0.309 + 0.099 







Aa = 


-2° 


2.303 + 0.017 


1.082 ± 0.015 


7.503 ± 0.077 0.863 ± 0.056 9.802 ± 0.045 


3.134 + 0.020 


0.317 + 0.104 




6.54 


Aa = 


-5° 


2.333 + 0.017 


1.087 + 0.015 


7.491 + 0.075 0.857 + 0.056 9.831 +0.045 


3.121 + 0.019 


0.389 + 0.104 




260.74 




FIG. 2: The one-dimensional posterior distributions of the tensor to scalar ratio r derived from the future Planck mock data 
without including the rotation angle. The fiducial values of r are (left panel) and 0.3 (right panel). The rotation angles 
considered in the fiducial model are: Aa = deg (black solid lines), Aa = —2 deg (red dashed lines), and Aa = —5 deg (blue 
dash-dotted lines). 



angle by a factor of 15, namely a(Aa) ~ 0.06 deg, which is consistent with our previous results [1,[T3- The non-zero 
rotation angle and the possible cosmological CPT violation can be verified by the future CMB measurements. 

However, if we neglect the effect of rotation angle at all (Aa = 0), the contribution of non-zero rotation angle on the 
BB power spectrum will be wrongly considered as the primordial tensor perturbations. For the fiducial model r = 0, 
the BB power spectrum due to the primordial tensor B-mode should vanish. But in the presence of non-zero fiducial 
rotation angle (Aa = —2, —5 deg), the BB power spectrum should be no n- vanishing. If we force the rotation angle to 
be zero in the analysis, the value of r will be enlarged to match the mock non-zero BB power spectrum. Using the Aa- 
free cases as the reference model, in Table [TTTI we can see that the minimal x 2 of the cases with Aa = are extremely 
larger, which implies that this kind of model does not fit the mock data well. In the left panel of FigJH we find that 
the larger the fiducial value of Aa is, the larger the obtained central value of r becomes. Consequently, this large value 
of r will enhance the CMB TT power spectrum at large scales. The contribution of primordial scalar perturbations 
on the TT power spectrum will be suppressed, which will change constraints of other cosmological parameters, for 
example, the baryon energy density Qt,h 2 becomes larger and the amplitude of primordial scalar perturbations A s 
becomes smaller. See the Table IIIII for details. Therefore, if we do not take the effect of rotation angle into account 



7 



0.0010 




Multipole / Multipole I 



+ 



r = 0.05 
r = 0.1 
r = 0.3 
r = 0.5 



100 1000 
Multipole / 



+ 



0.000 


I i i 


■?. \ /// 

•\ - /; 


\v 


0.001 


-•■}. i: 

■\ i! 




0.002 


'■■ \ i i 
! . v /' : 
" \ \y : 


r = 0.05 

r = 0.1 


0.003 




r = 0.3 

r = 0.5 ■ 


0.004 


1 




10 


100 1000 



Multipole / 



FIG. 3: Upper panels: The theoretical TB and EB power spectra for different values of e with r = 0.1. Below panels: The 
theoretical TB and EB power spectra for different values of r with e = — 1. We also show the theoretical TB power spectrum 
with the rotation angle Aq = —0.1 deg for comparison in the upper left panel (magenta dotted line). 



properly, the constraints of cosmological parameters derived from the future CMB data with high accuracy will be 
biased. We also find the similar situation for the case with r = 0.3 (right panel of Figj2]) . 



B. Gravitational Chern-Simons Term 



In this subsection, we study the effect of gCS term on the CMB power spectra and constrain the parameter e from 
the CMB polarization data. Note that, in our analyses we do not follow Ref.[12| which took the effects of eCS and 
gCS terms into account simultaneously, since the current constraint on e is still too weak to consider this kind of 
degeneracy. We always assume the rotation angle Aa = and neglect the systematic effects of CMB polarization 
data in the calculations. Therefore, we only have one more free parameter e G [—1,1], induced by the gCS term, and 
assume that it is scale- independent. 

Firstly, we show the effect of the gCS term on the CMB polarization power spectra. As we discuss before, a non-zero 
e only generates the non-zero TB and EB power spectra and leaves other four power spectra unchanged. In the upper 
panels of FigJ3] we plot the CMB TB and EB power spectra for different fiducial values of e with r = 0.1. We find that 
the non-zero e mainly affect the TB power spectrum on large scales {I < 200). There is a significantly enhancement 
on the amplitude of TB power spectrum on the largest scales (£ < 10), which is related to the CMB reionization 
information [2l[. We also notice that there are only two crossing points around £ ~ 15 and 70 on the TB power 
spectrum. For comparison, we show the TB power spectrum generated by a rotation angle Aa = —0.1 deg in the 
plot (magenta dotted line). As we know, the non-zero rotation angle will generate the TB power spectrum, based 
on the relation: C g TB = Cj E sin(2Aa) which implies that the shape of TB power spectrum is exactly the same with 
that of TE power spectrum. Therefore, the TB power spectrum with a non-zero Aa has many crossing points at high 
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FIG. 4: Marginalized one-dimensional and two-dimensional distributions (1, 2 a contours) of e and r. In the right below panel 
we also show the constraint of r without considering the non-zero gCS term, e = (red dashed line). 



multipoles (I > 200) and its signal mainly comes from the small scales. We could distinguish the TB power spectrum 
from the eCS term and that from the gCS term clearly on small scales. Moreover, in the plot we only show the TB 
power spectrum with Act = —0.1 deg which is much smaller than the best fit value from WMAP7+B03+BICEP data, 
Aa ~ — 2 deg, because the signal of TB power spectrum from the eCS term is much higher than that from gCS term, 
especially on small scales. Thus, if we find very large signal on small scales of TB power spectrum from the future 
observation, it only could be the signature of a non-zero rotation angle. 

We find the similar situation on the EB power spectrum. From the upper right panel of Figj3l we can see that the 
effect of the EB power spectrum, induced by the gCS term, is mainly on large scales and becomes to zero at high 
multipoles. There are many crossing points on the EB power spectrum, while the EB power spectrum, induced by 
the eCS term, does not have the crossing point. Based on the relation: C e EB ~ (Cf E — Cf B ) sin(4Aa), the EB power 
spectrum can not cross the zero line, since the EE power spectrum is always positive and larger than the BB power 
spectrum. Here, we do not show the curve in the plot for comparison, since the amplitude of EB power spectrum 
produced by the eCS term is much larger than that by the gCS term. Thus, the EB power spectrum can also be used 
to distinguish the eCS and gCS models, if the future CMB polarization data are accurate enough. 

In the below two panels of Figj3l we show the dependence of r on the TB and EB power spectra with e = — 1. One 
can see that the amplitudes of TB/EB power spectra are highly related to by the fiducial value of r, and their shapes 
remain unchanged. Therefore, the significance of detection of non-zero e will be directly determined by the primordial 
tensor perturbations r. At present, CMB observations do not detect the signature of non-zero r, which means that 
results on r are consistent with zero. Consequently, we can not constrain the parameter e very well from the current 
CMB polarization data, due to the very small primordial tensor perturbations. In FigUl we show the constraints 
on e and r from the WMAP7+B03+BICEP+QUaD data combination after marginalized over other cosmological 
parameters. From the one-dimensional distribution of e, as well as the fiat two-dimensional contours, apparently the 
present CMB data can not give good constraints on e. There is a peak around e = — 1, which implies that a right-hand 
polarized gravitational wave might be slightly favored. We also show the one-dimensional distribution of r including 
and not including the free parameter e. Due to the degeneracy between e and r, the constraint on r is slightly relaxed 
when a free e is included in the calculation. 

Thus, we use the same fiducial model in the last subsection to simulate the future Planck polarization data to 
constrain e. For the tensor to scalar ratio r, we have two fiducial models: r — 0.1 and r — 0.3. For each model, we 
consider six fiducial values of e with e — 0, —0.2, —0.4, —0.6, —0.8, —1. In Table HVl and upper panels of FigJS] we 
show the constraints on e from the future CMB mock data with different fiducial models. If we choose the fiducial 
value r = 0.1, even the Planck mock data can not give good constraint on e (the upper left panel). The 95% confidence 
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FIG. 5: The one-dimensional posterior distributions of e derived from the future Planck (upper panels) and CMBPol (below 
panels) measurements with different fiducial models of r. For each case, we have six fiducial values of e: e = (black), —0.2 
(red), —0.4 (blue), —0.6 (magenta), —0.8 (green), and —1 (yellow). 



TABLE IV: Constraints on e (95% C.L.) from future CMB mock data. 



Fiducial e 


Planck, r = 0.1 


Planck, r = 0.3 


CMBPol, r = 0.1 


CMBPol, r = 0.3 


e = 


[-1,1] 


[-0.69,0.66] 


[-0.71,0.70] 


[-0.43,0.42] 


e = -0.2 


[-1,1] 


[-1,0.57] 


[-1,0.54] 


[-0.61,0.23] 


e = -0.4 


[-1,1] 


[-1,0.32] 


[-1,0.37] 


[-0.80,0.05] 


e = -0.6 


[-1,1] 


[-1,0.22] 


[-1,0.25] 


[-1,-0.12] 


e = -0.8 


[-1,1] 


[-1,0.02] 


[-1,0.03] 


[-1,-0.36] 


e = -1.0 


[-1,1] 


[-1,-0.12] 


[-1,-0.08] 


[-1,-0.48] 



levels of e are always [—1,1] for six fiducial values of e. As we discuss before, the larger the value of r is, the higher 
the amplitudes of CMB TB/EB power spectra are (the below panels of Fig[3]). When the fiducial value of r becomes 
larger (r = 0.3), the constraints on e become better. The standard deviations of e have been shrunk significantly. But 
the constraints are still not conclusive. Obviously, the future Planck measurement will be very difficult to constrain e, 
if the primordial tensor perturbations are small, r < 0.1. Therefore, we also simulate the future CMBPol polarization 
data with higher accuracy and re-do the calculations, see Table IIIII and below panels of FigJSJ In this case, even 
with r = 0.1 the CMBPol data could also give good constraints on e. Finally, we obtain the standard deviation of e: 
a(e) ~ 0.1 with the fiducial value r = 0.3. The future CMBPol experiment could constrain e better. 

IV. CONCLUSIONS AND DISCUSSIONS 

Probing the violation of fundamental symmetries is an important way to search for the new physics beyond the 
standard model. In this paper we present constraints on the eCS and gCS models using the latest CMB polarization 
data, as well as the future simulated mock data. 
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In the constraints of rotation angle, induced by the eCS model, we extend our previous works by including the 
systematic errors of CMB polarization measurements. We consider two rotation angles, Aa and /3, in order to take into 
account the real rotation signal and a systematic error for each CMB polarization measurement, and impose priors 
on the systematic errors. Adding the systematic effects of CMB polarization data, we do not find significant change 
on the constraints of Aa, except the error bars become slightly larger. WMAP7+B03+BICEP data combination still 
favors a non-zero rotation angle about 2.2 a confidence level, namely Aa = — 2.28± 1.02 (deg), no matter whether the 
systematic effects are included. We still find a ~ 2 a tension between QUaD and WMAP7+B03+BICEP observations. 
When combining all CMB polarization data together, we obtain the tight constraint on the rotation angle at 95% 
confidence level: —1.34 < Aa < 0.82 (deg), which is consistent with a CPT-conserving Universe. 

Since the current constraints on the rotation angle are not conclusive, we simulate the future Planck polarization 
data. We find that the future CMB data could significantly improve the constraint of Aa by a factor of 15 (cr(Aa) ~ 
0.06 deg). Because a non-zero rotation angle will generate the BB power spectrum, we need to take this effect into 
account properly in the future data analysis. Otherwise the constraints of cosmological parameters will be apparently 
biased. 

We also give the constraint on the parameter e, which is induced by the gCS term and denotes the difference between 
the right- and left-handed polarized components. Since the effects of non-zero e on the TB and EB power spectra 
are very small, the parameter e is almost unconstrained from the current polarization data. Therefore, we simulate 
the future Planck and CMBPol data to constrain e further. We find that the future Planck data is very difficult to 
constrain e, if the primordial tensor perturbations are small, r < 0.1. Using the future CMBPol data, the constraint 
of e can be improved to cr(e) ~ 0.1, if the fiducial value r = 0.3. We could use the future CMBPol data to verify the 
gCS model. 
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